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Abstract. - Multi-electron giant dipole resonances of atoms in crossed electric and magnetic 
fields are investigated. Stationary configurations corresponding to a highly symmetric arrange- 
ment of the electrons on a decentered circle are derived, and a normal-mode stability analysis 
is performed. A classification of the various modes, which are dominated either by the mag- 
netic or Coulomb interactions, is provided. A six-dimensional wave-packet dynamical study, 
based on the MCTDH approach, is accomplished for the two-electron resonances, yielding in 
particular lifetimes of more than 0.1 fj,s for strong electric fields. 



Introduction. - Atoms exposed to strong external fields have proven to be a persistent 
source of intriguing phenomena with major impact on a variety of other fields such as the 
quantum dynamics of finite systems or Laser spectroscopy of Rydberg atoms [1-3]. Focusing 
on crossed electric and magnetic fields, new configurations of the hydrogen atom, the so-called 
giant dipole states, were discovered in the nineties [4-6]. These states represent a new form 
of microscopic matter where the electron and proton are separated by large distances due to 
the presence of an electric field: There exists an outer potential well whose bound atomic 
states possess a huge electric dipole moment. The atomic wave function then resembles an 
oscillating barbell in contrast to the usual shell-like structure of the electronic states in field- 
free space. In refs. [7] the giant dipole states of Positronium have been investigated, thereby 
arriving at the conclusion that the decentered positron-electron configuration is a quasistable 
state, i.e., the matter-antimatter system could be prevented from annihilation for very large 
time scales up to several years. Recently, dipolar matter has also become of major interest in 
the context of ultracold atomic and molecular physics, such as dipolar quantum gases [8] or 
ultracold molecular Rydberg states [9]. 

In the present investigation we consider an N— electron atom (total mass M) where the 
electrons (m) and the nucleus, interacting via the Coulomb potential V, are subject to crossed 
electric and magnetic fields. The translational invariance yielding the conservation of the 
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total canonical momentum in the absence of the fields is now lost, which manifests itself in 
the appearance of the coordinate-dependent vector potential. For homogeneous fields the 
conserved (total) pseudomomentum K = Y^iLo^i with = p.; — e^A^ + e^B x r$, whose 
components commute due to the neutrality of the atom, is a conserved quantity associated 
with the center-of-mass (CM) motion of the atom [10-13]. It can be exploited to perform a so- 
called pseudoseparation of the CM motion. The latter was originally done [10-13] for a fixed 
gauge. It is only recently that gauge-independent pseudoseparations for one-electron [6] and 
many-electron systems [14] have been performed. The latter provided a generalized potential 
for the electronic motion, possessing major impact on the prediction of new structures and 
effects in crossed fields such as the decentered giant dipole states [6,14]. Our starting point is 
the pseudo-separated Hamiltonian 



where the first two terms and the last three terms represent the gauge-dependent kinetic 
energy of the electrons relative to the nucleus (r-i = fj({pj, A(r 3 )}), i,j = 1,...,N) and 
the generalized gauge-independent potential, respectively. The potential contains, besides the 
Coulomb interaction V, Stark and diamagnetic interaction terms that are responsible for the 
existence of an outer well and bound giant dipole states for sufficiently strong (external and / or 
motional) electric fields in the one-electron case. 

Little is known, however, about multi-electron giant dipole resonances (GDR) which are 
the subject of the present investigaton: In Ref. [14] a naturally limited analytical investigation 
of the decentered two-electron configuration was performed. It provided evidence, but no final 
conclusion, on the existence of highly symmetric decentered resonances. In the following we 
derive the electronic configurations corresponding to the decentered TV— electron giant dipole 
resonances. We classify the modes and investigate their stability. Moreover, a numerically 
exact six-dimensional wave-packet dynamical investigation of two-electron resonances will be 
performed, supplying valuable information on lower bounds of the lifetimes of the GDR. 

Stationary decentered configurations. - With the gauge-independent generalized poten- 
tial at hand, let us search for decentered stationary points, being candidates for configurations 
of GDR. Inspecting the generalized potential (see Eq.Q, it is natural to introduce the elec- 
tronic center of mass (ECM) R = -h YliLi Y i as a new coordinate. In addition, we choose 
N ~1 vectors relative to the ECM, i.e., we decompose =: R+Sj, i = 1, . . . , N — 1. Defining 
Q — (R, si, . . . , sat_i) t we seek the stationary configuration 



Employing the coordinate frame £ = B/£>; 77 = TZ'/K'; £ = 77 x £/|tj x £| (where K' = 
K + AfE x B/_B 2 , and we assume Z(K',B) = 90°), one can show that there exist solutions 
that fulfill |rf 5 1 = |R + s| 0) | =: r Vi such that the ECM vector R = X£ + Yr] + ZC, is aligned 
along the £ axis, Rq = (Xq,0,0) t , where the decentering coordinate Xq satisfies 




(1) 




(2) 



Pr(Xo) 



NeK'B (NeB) 2 
M + M 



X + NZe 2 -^- = 0. 



(3) 



S. ZOLLNER et al: MULTI-ELECTRON GIANT DIPOLE RESONANCES IN CROSSED FIELDS 



3 



1500-, 



100IK 




Figure 1 - A giant-dipole configuration for N = 4 electrons (B = 1CT 4 , K = 2K CI ). The ECM 
R = (X, 0, 0) is decentered, with the relative vectors Si confined to a circle as indicated. Also shown 
is the electronic vector relative to the nucleus = R + s^. Typical motions for the cyclotron and 
Coulomb mode are indicated. 



Moreover, the relative coordinates are arranged on a circle in the orthogonal complement, 
Sj = s (OjCOS^^sin^) 71 . For symmetry reasons, we demand that all relative coordinates be 
distributed uniformly on that circle 

<f> k = <5> N k + A, $ N = 2tt//V. (4) 

This procedure determines the decentered configuration only up to a global rotation by an 
angle A g [0, 2ir). However, it allows us to fix the ratios of \X \ and s, r 



/ 47V 

r = s? _ m — =: a N s =: a N \X \ (5) 
V Zjfc=i (1/sm— ) 

and thus to solve Eq. ©. The solutions are 

where 6 = arccos[2(A' cr / ' K'f - 1], provided that K' > K CI = |^v/^- For any B there is 
a critical value K CI for K' at which the decentering sets in. Experimentally, the value of K' 
can be controlled via the strength of the external electric field. 

To conclude, the stationary electronic vectors possess the orthogonal decomposition = 
X£ + s(0, cos^i, sin0i) T with the common decentering along £. They are confined to a highly 
symmetric circular configuration perpendicular to £. This circular configuration is determined 
up to an overall rotation A. A generic setup is sketched in Figure The extremal solution 
exists only if the effective pseudomomentum K' exceeds some critical value given above. With 
increasing K' the decentering of the ECM becomes more pronounced. 

Normal-mode analysis of N -electron giant dipole states. - We expect the decentered 
stationary configurations to be promising candidates for GDR with certain lifetimes. Pre- 
ceding a numerical study, we first seek to obtain some insight into the local stability of the 
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extremal configurations. According to Ehrenfest's theorem, the expectation values (Q)(t) of a 
harmonic system obey the corresponding classical equations of motion. This suggests treating 
the problem within a normal-mode analysis. 

The equations of motion are obtained in terms of the displacements v(t) :— Q(t) — Qo, 

v = u> ■ v + A ■ v. (6) 

The antisymmetric 3N x 3N cyclotron matrix u> contains the Lorentz force and the harmonic 
matrix A is built up essentially from the Hessian of the generalized potential. The solution of 
this system of differential equations is given by the span 

6^ 

v(t) = $> P e 7pt )cp (v P S C 3N ; Jp , c p e C), (7) 
P =\ 

fulfilling the quadratic eigenvalue equation 

( 7p 2 / - 7p w - A)v p = 0, (8) 

/ being the identity. Our stability analysis amounts to finding the complex eigenvalues 7 P =: 
T p + itt p (whose imaginary parts are frequencies of a vibration about a stable point, and 
whose real part corresponds to an instability), and the eigenvectors v p . The above quadratic 
eigenvalue problem is solved by reducing it to the standard linear eigenvalue problem 

Every eigenpair (7, v) of © with Inry 7^ has a twin pair (7*,?;*), Thus we treat the 6./V 
modes as effectively 3N modes. 

If one studies the behavior of the modes {^p(K] B)} p for different electron numbers N 
in the vertical configuration A = —ir/2, one observes that even though the patterns become 
increasingly rich and involved with larger N, there is a clear distinction between the character 
of the modes regarding their behavior as a function of K (and B) , their order of magnitude 
and, more generally, their location in the complex plane. Based on our analysis, including the 
associated eigenvectors, we arrive at the following classification of the 37V eigenmodes: 

• N so-called cyclotron modes corresponding to the cyclotron motion of the effective electronic 
particles. Their values are exclusively of the order of the cyclotron frequencies, 7 = id ~ i^^-, 
i.e., almost independent of K. 

• N Coulomb modes corresponding to the inter-Coulombic motion, with dominating contri- 
butions from the harmonic matrix A. They fall off quickly with K, since the decentering 
Xq(K) increases, and the Coulomb interaction weakens. As opposed to the cyclotron modes, 
the Coulomb motion takes places predominantly parallel to B. 

• 1 CM mode, roughly reflecting the cyclotron motion of the CM, Q ~ nnr- 

• 1 zero mode (7 P = 0) stemming from the rotational invariance of the saddle point with 
respect to the circular configuration (see Eq. 0} 

• N — 2 modes termed decay modes in recognition of the fact that they are predominantly 
real. They are neither directly related to the cyclotron motion nor to the spectrum of A, and 
their slope is twice as steep as that of the Coulomb modes. Their absence for the case N = 2 
is essentially why the two-electron system is locally stable. 

A generic example for the spectrum is given in Fig. for the case of N = 4 electrons. 
There are a few subtleties that go beyond the categorization suggested above. Without going 
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Figure 2 - Eigenmodes {7 P } P for N — 4 electrons as a function of the electric field E = BK/M, shown 
for the range K/K CT € [1, 10] (B = 1CP 4 ). Imaginary parts appear as solid lines, while points (+) are 
used for the real parts. - The top horizontal line represents four almost degenerate cyclotron modes. 
Below, the four Coulomb modes fall off quickly and intersect the CM mode (the nearly horizontal 
line about four orders below the cyclotron modes). There are two residual decay modes, whose slope 
is twice that of the Coulomb modes. 



into the details, we mention that there are certain interactions among different mode types. 
Their principal causes are crossings between the CM mode and the decay modes (resulting 
in some striking deformations of the usual line pattern), and avoided crossings of the CM 
mode with at least some of the Coulomb modes. This may be taken as a hint for the different 
symmetry relations among the Coulombic and the decay modes. 

Let us remark on the influence of the global rotations A. While the stationary character 
is not affected by a common rotation of the relative coordinates Sj by an angle A on the 
circle, the dynamics differs. In order to see this dependence, we inspect the modes as a 
function of the rotation angle (7 P (A)}. Apart from an obvious symmetry — note that rotating 
by <J>at = 2n/N gives an indistinguishable setup — we found that the Coulomb modes and 
the decay modes show a pronounced periodic change. For the case N = 2, we encounter an 
exceptional behavior whenever A comes close to mod n. One of the Coulomb modes then 
tends to zero, along the way turning real. In this respect, the local stability indicated above 
is applicable only outside the singular horizontal configuration A = 0. 

Wave-packet dynamical study. - We now turn to a numerical study of the two-electron 
system (N = 2). We emphasize that a six-dimensional resonance study is at the frontier of 
what is currently possible and requires a careful choice of the computational approach. This 
applies especially in view of the fact that our system is governed by dramatically different time 
scales (see below). Therefore we adopted the Multi-Configuration Time-Dependent Hartree 
(mctdh) method [15,16], a wave-packet propagation scheme known for its outstanding effi- 
ciency in high-dimensional settings. Its basic idea is to solve the time-dependent Schrodinger 
equation by expanding the wave function in a moderately sized time-dependent basis related 
to Hartree products. We stress that this approach is designed for distinguishable particles. 
Applying it to a fermionic system like ours finds its sole justification in the fact that the spatial 
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separation between the electrons is so large that they are virtually distinguishable. 

Let us first point out how the computational method can be applied to our problem. There 
are essentially two conflicting types of motion and corresponding scales: perpendicular to the 
magnetic field, we have the magnetic length Rb ~ 10 2 a.u. (for the strong laboratory field 
strength B = 10~ 4 a.u. we consider). For the Coulomb motion parallel to B — expected to 
take place approximately in a harmonic potential — we can use the usual oscillator lengths 
zq = 1/ v / tow7 as an estimate. These are roughly on the order of 10 3 a.u. but increase with 

K (see above). Analogously, the anticipated time scales are T B = ^ B ~° 1.5 ps for the 
cyclotronic motion and 1.5ns — 15ns for the Coulomb motion. 

Within our numerical study we focus on the investigation of the stability of the giant dipole 
states by using the propagation of harmonic-oscillator wave packets initially localized at the 
extremal position. Relaxation techniques have been applied in order to improve the intially 
chosen wave packet. Moreover, following the evolution of wave packets with an initial displace- 
ment from the extremum, we tested the robustness of the resonances. For comparison, we also 
examine the decentered eigenvectors. The chosen parameter set is K/K cr 6 {1.1, 2.0, 10.0} 
(we drop the prime here and in the following), which accounts for the cases of just above 
threshold, the medium range and the very large K regime. For simplicity, we first focus on 
the vertical configuration A = —ir/2. The propagation times are chosen in the regime of 
50 — 100 ns; this includes many periods of the Coulomb modes and some 10,000 periods of 
the rapid cyclotron motion ! 

Very close to the critical point, e.g. K = 1.1K CX , the observed motion is unstable for 
several degrees of freedom, specifically for those belonging to the relative vector sj_ = (x, y) T 
perpendicular to the magnetic field, on a time scale of 10 3 ps. Other modes are affected too 
via couplings. This is a general fact, but for K = l.lA' cr it is very pronounced, whereas it 
is suppressed for the time evolution in case of larger values of K due to the weaker Coulomb 
modes (see Fig. EJ, 

For K being twice its critical value, K = 2K cr , the vertical configuration is virtually 
stable on a scale of T ~ 10 4 ps. As opposed to the case K/K CT = 1.1, the instability is 
almost exclusively due to the relative motion in y. However, the response of the system 
upon displacing the initial wave packet in Z and z by 2, 000 a.u. unveiled that x is rendered 
rather unstable. Hence the resonance is expected to be less robust, although its lifetime is 
not significantly reduced altogether. To illustrate the resonance character of the system, Fig. 
El gives an impression of the excitation spectrum, obtained via Fourier transformation of the 
auto-correlation function c(t) :— (^o I e~ tHt ^o). The 'displaced state' with (z)q = (Z)q = 
2. 000 a.u. produces a rather interesting excitation spectrum. The equidistant spacing of the 
peaks can be interpreted as a signature of harmonicity for both excited degrees of freedom 
(Z, z). 

For K = 10K cr , we find the system to be practically stable on the time scale T = 10 5 ps. 
The initial wave packet experiences only tiny deformations with respect to all degrees of 
freedom but z, showing some oscillatory behavior. However, those are still marginal compared 
to the spatial extension of the decentered state. As a result of our wave-packet dynamical 
study, we can conclude that the lifetimes of the GDR for sufficiently large K are beyond 0.1/xs. 

To complete the discussion, let us touch on the effects of different circular configurations A 
of the GDR. As examples, we investigated both the supposedly unstable 'horizontal configura- 
tion' (A = 0) and a 'diagonal configuration' (A = — tt/4), corresponding to different settings 
of the extremal relative coordinates s^ ^ . To sum up our findings, the horizontal configuration 
indeed adds an instability, which is discernible even for very high K, if less distinct. The 
diagonal configuration was partly unstable on a timescale comparable to that of the vertical 
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Figure 3 - Spectrum for K/K cr = 2 for the excited initial state (displaced in Z, z). 
case, T ~ 10 4 ps. 

Conclusion and outlook. - Our investigation shows that atoms in crossed fields exhibit 
multi-electron giant dipole states with extraordinary lifetimes. These resonances constitute 
highly symmetric and exotic states of matter, where the elecrons are strongly correlated and 
can, for laboratory field strengths, be separated from the nucleus (or positively charged core) 
by many thousand Bohr radii. An experimental preparation of the giant dipole states might 
employ the scheme suggested for single electrons [17]. The latter is based on the preparation 
of Rydberg states via Laser excitation, followed by a sequence of electric-field switches that 
carry the excited electrons to the decentered configuration. Preparing different initial Rydberg 
states, like circular ones, enhances the variety of accessible decentered states. 
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